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Abstract 

We study the decays of top squarks (fi,2) and bottom squarks (61^2) in the Min- 
imal Supersymmetric Standard Model (MSSM) with complex parameters At,Aii,fi 
and Ml. We show that including the corresponding phases substantially affects the 
branching ratios of ti^2 and 61^2 decays in a large domain of the MSSM parameter 
space. We find that the branching ratios can easily change by a factor of 2 and 
more when varying the phases. This could have an important impact on the search 
for ti^2 and 61^2 and the determination of the MSSM parameters at future colliders. 



1 Introduction 



The experimental studies of supersymmetric (SUSY) particles will play an important 
role at future colliders. Studying the properties of the 3rd generation sfermions will be 
particularly interesting because of the effects of the large Yukawa couplings. The lighter 
of their mass eigenstates may be the lightest charged SUSY particles and they could be 
investigated at an e~^e~ linear collider p. Moreover, they could also be copiously produced 
in the decays of heavier SUSY particles. Several phenomenological studies on SUSY 
particle searches have been performed in the Minimal Supersymmetric Standard Model 
(MSSM) 12] with real SUSY parameters. Analyses of the decays of the 3rd generation 
sfermions ti^2, &i,2; ^1,2 and u^- in the MSSM with real parameters have been made in 
Refs. i| , and phenomenological studies of production and decays of the 3rd generation 
sfermions at future e~^e~ linear colliders in Ref. [3]. 

In general, however, some of the SUSY parameters may be complex, in particular 
the higgsino mass parameter /i, the gaugino mass parameters Mi 2,3 and the trilinear 
scalar coupling parameters ^/ of the sfermions /. The SU(2) gaugino mass parameter 
M2 can be chosen real after an appropriate redefinition of the fields. The experimental 
upper bounds on the electric dipole moments (EDMs) of electron, neutron and the ^^^Hg 
and atoms may pose severe restrictions on the phases of the SUSY parameters, 

though they are model dependent. In the constrained MSSM the phase of fj, turns out 
to be restricted to the range \ip^\ < 0.1 if all SUSY masses are in the TeV range 0. 
On the other hand, in more general models, such as those with lepton-flavour violation, 
no restriction on the phase ip^ from the electron EDM at one- loop level is obtained [?]. 
Restrictions arising due to two-loop contributions to the EDMs are less severe jS]. 

Therefore, in a complete phenomenological analysis of production and decays of the 
SUSY particles one has to take into account that Af,fi and Mi can be complex. The most 
direct and unambiguous way to determine the imaginary parts of the complex SUSY pa- 
rameters could be done by measuring relevant CP-violating observables. For example, for 
the case of sfermion decays CP- violating jUl 02] and T- violating ^I] observables have been 
proposed. On the other hand, also the CP-conserving observables depend on the phases 
of the complex parameters, because in general the mass-eigenvalues and the couplings 
of the SUSY particles involved are functions of the underlying complex parameters. For 
example, the decay branching ratios of the Higgs bosons depend strongly on the complex 
phases of the i and b sectors [121 HSl El , while those of the staus fi 2 and r-sneutrino 
can be quite sensitive to the phases of the stau and gaugino-higgsino sectors [T3| . Also the 
Yukawa couplings of the third generation sfermions are sensitive to the SUSY phases at 
one-loop level JH] ■ Furthermore, explicit CP violation in the Higgs sector can be induced 
by i and b loops if the parameters At, Af, and /i are complex [T21 [13 [TH]. It is found 
[121 [m dZl [Hj that these loop effects could significantly influence the phenomenology of 
the Higgs boson sector. 

In this article we study the effects of the phases of A^, Af,, and Mi on the decay 



2 



branching ratios of the ti^2 and 61,2 with qi (^2) being the hghter (heavier) squark. We 
take into account the exphcit CP violation in the Higgs sector. We will show that the 
influence of the phases can be quite strong in a large domain of the MSSM parameter 
space. This could have an important impact on the search for ti_2 and 61 2 and on the 
determination of the MSSM parameters at future colliders. 

In section 2 we discuss the SUSY CP phase dependences of masses, mixings and 
couplings. Section 3 contains our numerical investigation on the CP phase dependence of 
the branching ratios of ii^2 and 61^2 decays. In section 4 we present our conclusions. 



2 SUSY CP Phase Dependences of Masses, Mixings 
and Couplings 

In the MSSM the squark sector is specified by the mass matrix in the basis {qi, qpt) with 
g = t or 6 |2ni 

^2 ^ ( < aim, \ 



agmq mt^ 



with 



= M? +m| cos 2/? (J|^ -eg sin^^VK) (2) 
= M0 + m| cos 2(3 sin^ 6^ + m^, (3) 

_ f {A- fi* cot P)mt {q = i) 
" " \ (A - /i* tan/5) m, {q = b) ^ ' 

= \aqmq\ e"^^ (-vr < (pq < vr). (5) 

Here /f is the third component of the weak isospin and Cq the electric charge of the quark 
q. Mq fj and Af h are soft SUSY-breaking parameters, /x is the higgsino mass parameter, 
and tan/5 = V2/V1 with vi (^2) being the vacuum expectation value of the Higgs field 
{H2). As the relative phase ^ between vi and V2 is irrelevant in our analysis, we adopt the 
^ = scheme jTTj. We take Aq (g = t,b) and /i as complex parameters: Aq = \Aq\ e^'^^i 
and /i = |/i| e*'^'" with — vr < ^Aq,^l < tt. Diagonalizing the matrix one gets the mass 
eigenstates qi and q2 

pg( \ ^( e'^^cosOq sinOq \ ( Ql \ ..x 
[qr J [ -sinOq e-'^^cosOq J [ qn J 

with the masses m^^ and mq^ (m^^ < m^j), and the mixing angle 9g 

^1,2 = + ^Ir T \/ ("^i - + 4|agmgP), (7) 

Oq = te.n-\\aqmq\/{ml-mlJ) {-7r/2 < Oq < 0). (8) 




3 



The Ql — Qr mixing is large if |m?^ — m?^| < |agmg|, which may be the case in the t sector 
due to the large rut and in the b sector for large tan /3 and From Eqs. (jH),© and (jHl) 
we see that m?^ ^ and 9q depend on the phases only through a term cos{(pAq + V'^i)- This 
phase dependence of the i (6) sector is strongest if \At\ ^ cot (3 {\Ah\ ^ tan/3). 

The properties of charginos xf {i = 1,2; m-± < fn^±) and neutralinos x° (j = 1, 4; 
m^o < ... < wi^o) are determined by the parameters M2, Mi, fi and tan /5. We assume that 
the gluino mass rrig is real. We write the U(l) gaugino mass Mi as Mi = iMile*"^^ (— vr < 
(pi < 7i). Inspired by the gaugino mass unification we take |Mi| = (5/3) tan^ ^vkM2 
and nig = {asirrig) / a2)M2. In the MSSM Higgs sector with exphcit CP violation the 
mass-eigenvalues and couplings of the neutral and charged Higgs bosons, H^, H2, 
{mffo < mjjo < niffo) and if^, including Yukawa and QCD radiative corrections, are 
fixed by 112^+, tan P, fi, nit, rrib, Mq, M^^, Mfj, At, Ah, \Mi\, M2, and rrig [T7]. The neutral 
Higgs mass eigenstates ifj* , H2 and are mixtures of CP-even and CP-odd states (0i,2 
and a) due to the explicit CP violation in the Higgs sector. For the radiatively corrected 
masses and mixings of the Higgs bosons we use the formulae of Ref . ^Tj . 

Here we list possible important decay modes of ti,2 and bi^2'- 

h ^ t~g,tx'^,bxl,b,W+,hH+ (9) 

k ^ t~g,tx^,bxl,hZ\h,2W+,hH',,\2H+ (10) 

bi ^ b~g,bx1,tx] ,iiW- ,hH- (11) 

62 ^ bg,bxltxJ,hZ\h,2W-,hHlh,2H-. (12) 

The decays into a gauge or Higgs boson in ©-(1121) are possible in case the mass splitting 
between the squarks is sufficiently large ^ |4j . The explicit expressions of the widths of 
the decays ©-(^21) in case of real SUSY parameters are given in [21j. Those for complex 
parameters can be obtained by using the corresponding masses and couplings (mixings) 
from Ref s. [T71 |2ni 1211 and will be presented elsewhere |221- 

The phase dependence of the widths stems from that of the involved mass-eigenvalues, 
mixings and couplings among the interaction-eigenfields. Here we summarize the most 
important features of the phase dependences. 

(I) qi (g = i, b) sectors: 

(a) The mass-eigenvalues rrig-^^ ^ are sensitive to the phases {ipAq,^^) via cos(v9^^ -|- 
ip^) if and only if la^mgl ~ ("^^^ +"^^fl)/2 o-n-d \Ag\ ~ \fJ'\Cg (with Ct = cot P 
and Cb = tan P). 

(b) The g- mixing angle 9g (given by tan26'q = 2|agmg|/(m?^ — "m-?^) ) is sensitive 
to (99^,, V^A*) via cos(v9a, + V^^) if and only if 2\agmg\ > |m?^ - m?^| and |v4q| ~ 

\fl\Cg. 
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(c) The g-mixing phase ipq in Eq.® is sensitive to (v^a,, V'^t), V^a, and y?^ if |y4g| ~ 
|/i|Cg, \Aq\ ^ \fJ'\Cq and \Aq\ <^ |yu|C*g, respectively. For large squark mixing 
the term cx: sin 26^^ cos yjg can result in a large phase dependence of the decay 
widths (see Eq.®). 

(II) Xi and xf sectors: 

(a) m^o {i = 1,...,4) and the mixing matrix are sensitive [insensitive] to the 
phases {ipi,ip^) for small [large] tan/?. 

(b) m-± and the x^-mixing matrices are sensitive [insensitive] to <y9^ for small 
[large] tan/3. 

(III) sector: 

This sector is independent of the phases. 

(IV) sector: 

(a) m^o {i = 1,2,3) are sensitive [insensitive] to the phase sums '^a^ + V^ai ^^t 
small [large] tan/5 [T7j . 

(b) In general the if^—mixing matrix (a real orthogonal 3x3 matrix Oij) is sensitive 
to ipAt b + ffi for any tan P |T7j. 

(V) The couplings among the interaction-eigenfields: 

(a) For the decays into fermions and gauge bosons in Eqs.(P)- (fT^ . they are gauge 
couplings and/or Yukawa couplings (ht^), and are independent of the phases 
at tree-level. 

(b) For the decays into Higgs bosons in Eqs.(P|)- (fT^ . the qL-QR {ql-Ql, Qr-Qr) 
couplings are dependent on (independent of) the phases ipAt, '^a^ and y^^; 



with 



C(h{tRH-) 


~ sin P ht{A1 cot P + fi) 


(13) 


C{AbRH^) 


~ cos P hb{Al tan P + fi) 


(14) 




~ htfi 


(15) 


C(fe02) 


~ htA; 


(16) 


C{titRa) 


~ sin P ht{ a; cot P + fi) 


(17) 




~ hAi 


(18) 




~ hbfJ^ 


(19) 


Cib{bRa) 


~ cos P hb{Al tan P + fi) 


(20) 


ht = 


gmt/{V2mw sin/5) 


(21) 


hb = 


gmb/{V2mw cos/?). 


(22) 
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Here 0i = Oijif° (i = 1,2) and a = OsjH^ are the CP-even and CP-odd 
neutral Higgs bosons, respectively [T7] . 

According to items (I)-(V) we expect that the widths (and hence the branching ratios) of 
the decays ©-(US) are sensitive to the phases (v^a*, V^a^, V^^i, V^i) in a large region of the 
MSSM parameter space. 

3 Numerical Results 

Now we turn to the numerical analysis of the ti^2 and 61,2 decay branching ratios. We 
calculate the widths of all possibly important two-body decay modes of Eqs.Q-lO- 
Three-body decays are negligible in the parameter space under study. In order to improve 
the convergence of the perturbative expansion I2S] we calculate the tree-level widths 
by using the corresponding tree-level couplings defined in terms of "effective" MSSM 
running quark masses m^'^" (i.e. those defined in terms of the effective running Yukawa 
couplings /ij'^" oc nf^^). For the kinematics, e.g., for the phase space factor we use 
the on-shell masses obtained by using the on-shell (pole) quark masses Mt^. We take 
Mt = 175 GeV, Mi, = 5 GeV, m™" = 150 GeV, m);™ = 3 GeV, mz = 91.2 GeV, 
sin^ = 0.23, m^y = nizcosOw, Oi{mz) = 1/129, and as{mz) = 0.12 (with as{Q) = 
Vl-n j ((33 — 2nj) ln(Q^/A^^)), being the number of quark flavors). In order not to vary 
too many parameters we fix \At\ = \Ab\ = \A\ and M2 = 300 GeV, i.e. rrig = 820 GeV. 
In the following we will assume that rrig > m^^ so that the decays ti — * tg and hi —>■ bg 
are kinematically forbidden. In our numerical study we take tan/?, m^^, m^^, 'm,^^, |y4|, 
(fAt: V'Ab, ViJ.y V'l &nd 771^+ as input parameters, where mi_^^ and are the on-shell 
squark masses. Note that for a given set of the input parameters we obtain two solutions 
for (Mq, Mjj) corresponding to the two cases m^-^ > m^^ and m^-^ < m^^ from Eqs. 
(^-(^ and (jZj) with mt replaced by Mt. In the plots we impose the following conditions 
in order to respect experimental and theoretical constraints: 

(i) m-± > 103 GeV, m^o > 50 GeV, m^^ j,^ > 100 GeV, m^^ > m^^o, 
m^o > 105 GeV, 

(ii) I < 3 (M? + M? + ml), and lA^p < 3 (M? + M? + m?), where mj = {mjr+ + 
ml sin^ B^) i5-\ml and ml = (mj,, + ml sin^ 6^) cos^ /3 - \ ml, 

(iii) Ap{i-b) < 0.0012 [21 using the formula of fIF] . 

(iv) 2.0 ■ 10^^ < B{b s'-f) < 4.5 ■ 10~^ I^Hl assuming the Kobayashi-Maskawa mixing 
also for the squark sector. 

Condition (i) is imposed to satisfy the experimental mass bounds from LEP |j27j . (ii) is the 
approximate necessary condition for the tree- level vacuum stabihty |2H1- (hi) constrains 
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H and tan/3 (in the squark sector). For the calculation of the b ^ S'j width we use the 
formula of j^H] including the 0{as) corrections as given in pUj . 

As alreaday mentioned the experimental upper limits on the EDMs of electron, neutron, 
^^^Hg and ^^^Tl strongly constrain the SUSY CP phases jH]. One interesting possibility for 
evading these constraints is to invoke large masses (much above the TeV scale) for the first 
two generations of the sfermions keeping the third generation sfermions relatively 
light ( < 1 TeV). In such a scenario [cpi, ip^) and the CP phases of the third generation 
[ifiAt^ VAf,^ ^Ar) are practically unconstrained We adopt this scenario. Furthermore, 
we have checked that the electron and neutron EDM constraints at two-loop level jH] are 
fulfilled in the numerical examples studied in this article. 

In Fig.l we plot the contours of the branching ratios of the ti decays B(ti txi) and 
B(ti — ^ bxf) in the \A\ — plane for ipAt = and tt/2 with tan/? = 8, {m^^,m^^,mij = 
(400,700,200) GeV, ipA^ = V^i = 0, = vr, and mH+ = 600 GeV in the case m^^ > m^^. 
In the case m^-^ < m^^ we have obtained similar results. As expected, these branching 
ratios are sensitive to the phase ip^t in a sizable region of the \ A\ — plane. The difference 
between (pAt = (Figs, la and c) and (pAt = vr/2 (Figs, lb and d) can be explained by 
item (I). Especially the strong (pAt dependence of the decay width T{ti bxt) is a result 
of item (I)(c) Moreover, the t-mixing angle 9i is sensitive to cos{(pAt + V^^i) for large 
|y4t|(= |y4|) with \At\ ~ |Ai|/8, which is a consequence of item (I)(b). 

In Fig. 2 we plot the contours of the ti decay branching ratios B{ti — *• txi) and 
B{ii bxt) in the ipAt - plane for tan/3 = 8, [mi^.mi^, m^J = (400, 700,200) GeV, 
\A\ = 800 GeV, = 500 GeV, 9?Ab=V^i=0, and mH+ = 600 GeV in the case m^^ > m^^. 
In the case mz < mx we have obtained similar results. One sees that these branching 
ratios depend quite strongly on the CP phases ipAt and (p^, as follows from item (I). 

In Fig. 3 we show the (p^ dependence of the ti decay branching ratios B{ti txi) 
and B{ii bxt) for = and tc/2 with tan/3 = 5, {mi^.mi^, m^J = (400,700,200) 
GeV, \A\ = 800 GeV, = 500 GeV, ipA, = ^a, = 0, and mH+ = 60o'GeV in the case 
THf > rrif . For mi < mz we have obtained similar results. As can be seen in Fig. 3 
the branching ratios are very [somewhat] sensitive to ip^ [yji] for small tan/?=5. For this 
choice of parameters the masses and mixings of t, 6, x° and are sensitive to ip^, while 
only the masses and mixings of x° are sensitive to (pi , as can be seen from items (I) and 
(II). In general, according to items (I), (II) and (V), the ti decay branching ratios tend to 
be sensitive [insensitive] to (pAt [v^a^, V^/^, and ipi] for large tan/?( > 15) in case the width 
of the biH^ mode in Eq. Q (which can be sensitive to (v^At, V^aJ as seen from Eqs. (0 
IT^ I14p ) is relatively small. We have confirmed this. 

In Fig. 4 we show the (pAt dependence of the ^2 decay branching ratios for y?^ = and 
tt/2 with tan/3=8, (mi^,mi^ ,m^J = (400,700,200) GeV, \A\ = 800 GeV, = 500 GeV, 
(pA^=<Pi=0, and mH+ = 600 GeV in the case m^^ > m^^. The case m^^ < m^^ leads to 
similar results. We see that the ^2 decay branching ratios are very sensitive to ipAt and 
depend significantly on yj^. This can be expected from items (I), (II), (IV) and (V)(b). 
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From item (I)(c) [(IV) (b) and (V)(b)] we expect that the widths of the decays t2 bxi,2 
[tiH^ and in Eq. lfTIH) can be sensitive to (pAt ['^Att and y?^] for large tan/9( > 15) 

if they are kinematically allowed. We have confirmed this. 

For 61 2 decays we have obtained results similar to those for the ti^2 decays. Here we 
show just a few typical results for them. In Fig. 5 we show the (f^^ dependence of the bi [62] 
decay branching ratios for tan/5=30, {m^^, m^-^, m^J = (200,700,400) GeV [(175,500,350) 
GeV], (|A|, |/x|) = (800,700) GeV [(600,500) GeV], ipA,=^^^=^T, <^i=0, and mH+ = 180 GeV 
in the case m^^ > m^^. For m^^ < m^^ our results are similar. We find that the 61,2 decay 
branching ratios are very sensitive to (pA^ expected from items (IV) (b) and (V)(b). The 
main reason is that the decay widths for 61^2 ~^ tiH~ strongly depend on ^pA^ (and y?^) 
for large tan/3 (see Eqs. (ll3ll4p ). This explains also the tendency of the ipA^ dependence 
of the branching ratios for the other decays shown in Fig.5. For small tan/5 ~ 8 we 
expect that the 61,2 decay branching ratios can be somewhat sensitive to ipA^ ^ and ip^ (see 
items (I) and (V)(b)). Similarly, we expect that the decay widths of ti,2-f^~ can be fairly 
sensitive to ipAt ^ and ip^. We have confirmed this. 

4 Conclusions 

We have calculated the branching ratios of the two-body decays of ii^2 and 6i_2 and studied 
their CP phase dependence within the MSSM with complex parameters At, Ab, fi and Mi. 
We have shown that the effect of the SUSY CP phases on the branching ratios can be 
quite strong in a large domain of the MSSM parameter space. The ipAt dependence of the 
61^2 decay branching ratios is mainly due to the ipAt dependence of their couplings to the 
Higgs bosons. In the case of the ii^2 decays the branching ratios depend on ipAt also via 
the t-mixing phase ip^ ~ ipAt for |^t| ^ 1/^1/ tan/5, in addition to the (pAt dependence of 
the mixing angle % and their couplings to the Higgs bosons. Some of the branching ratios 
can change by a factor of 2 or more when varying the phases. This CP phase dependence 
of the branching ratios could have an important impact on the search for ti 2 and 61^2 and 
on the determination of the MSSM parameters at future colliders. 
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Figure Captions 



Figure 1: Contours of the ii decay branching ratios B{ii txi) (a,b) and B(ii — * bxt) 
(c,d) in the \A\ — \fj,\ plane for (pAt — (a,c) and 7r/2 (b,d) with tan/3 = 8, (mf^,mf2,m^J 
= (400,700,200) GeV, (^^j=(/9i=0, (Pfj,=7i, and m//+ = 600 GcV in the case mf^ > JTi't^- 
The blank areas are excluded by the conditions (i) to (iv) given in the text. 

Figure 2: Contours of the ti decay branching ratios B{ti txi) (a) and B(ti — * bxf) 
(b) in the v^a^V^m Plane for tan/3 = 8, {mi^,mi^,mij= (400,700,200) GeV, \A\ = 800 GeV, 
= 500 GeV, (^A(,=<^i=0, and mH+ = 600 GeV in the case m^^ > m^^. 

Figure 3: dependence of the ii decay branching ratios B(ti txi) B{ti — hxt) 
for = (solid lines) and7r/2 (dashed lines) with tan /3 = 5, {mi^,m^^,mi ) = (400,700,200) 
GeV, \A\ = 800 GeV, = 500 GeV, ^At=^A,=0, and mH+ = 600 GeV in the case 

Figure 4: (pAt dependence of the ti decay branching ratios for = (a) and 7r/2 
(b) with tan/3=8, (mj^,mj2,m^J= (400,700,200) GeV, \A\ = 800 GeV, |//| = 500 GeV, 

¥'yij,=¥'i=0, and mii+ — 600 GeV in the case mf^ > mf^. Note that the t\H.2,-i 
modes are kinematically forbidden here. 

Figure 5: pAt, dependence of the 6i (a) [62 (b)] decay branching ratios for tan/3=30, 
(mt^,mt2,m^J= (200,700,400) GeV [(175,500,350) GeV], (|A|, |//|)=(800,700) GeV [(600,500) 
GeV], ipAt—^ii—T^-i </'i=0, and mH+ = 180 GeV in the case m^^ > mi^. Only interest- 
ing modes are shown in Fig.b where we have mg^ ~ 570GeV and {mfjOjmijo, 111^0) ~ 
(114,156,158) GeV. 
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Fig.5 
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